First-order coherence versus entanglement in a nano-mechanical cavity 
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The coherence and correlation properties of effective bosonic modes of a nano-mechanical cavity composed 
of an oscillating mirror and containing an optical lattice of regularly trapped atoms are studied. The system 
is modelled as a three-mode system, two orthogonal polariton modes representing the coupled optical lattice 
and the cavity mode, and one mechanical mode representing the oscillating mirror. We examine separately the 
cases of two-mode and three-mode interactions which are distinguished by a suitable tuning of the mechanical 
mode to the polariton mode frequencies. In the two-mode case, we find that the occurrence of entanglement 
in the system is highly sensitive to the presence of the first-order coherence between the modes. In particular, 
the creation of the first-order coherence among the polariton and mechanical modes is achieved at the expense 
of entanglement between them. In the three-mode case, we show that no entanglement is created between the 
independent polariton modes if both modes are coupled to the mechanical mode by the parametric interaction. 
There is no entanglement between the polaritons even if the oscillating mirror is damped by a squeezed vacuum 
field. The interaction creates the first-order coherence between the polaritons and the degree of coherence can, 
in principle, be as large as unity. This demonstrates that the oscillating mirror can establish the first-order 
coherence between two independent thermal modes. A further analysis shows that two independent thermal 
modes can be made entangled in the system only when one of the modes is coupled to the intermediate mode 
by a parametric interaction and the other is coupled by a linear-mixing interaction. 

PACS numbers: 03.67.Mn, 42.50.Ar, 42.50.Pq, 42.50.Wk 



I. INTRODUCTION 

The generation of quantum effects in nano-mechanical cav- 
ities with movable mirrors has been the subject of a great 
interest in recent years [1, 2]. This interest stems from the 
possibility of the development of new practical techniques for 
cooling macroscopic objects to very low temperatures and for 
engineering of entangled states of macroscopic systems. With 
the recent progress in laser cooling techniques, fabrication of 
low-loss optical elements and high-Q mechanical resonators, 
it is now possible to prepare nanomechanical oscillators that 
can be controlled to a very hight precision and can even reach 
the quantum level of the oscillations [3]. In these systems, the 
vibrations of mechanical oscillators are induced by radiation 
pressure that creates a strong nonlinear coupling of the vibra- 
tional mode to radiation modes. Most of these studies have 
been done on examples provided by cavity optomechanical 
systems with linear or ring cavities [4]. It has been demon- 
strated that a large radiation pressure can be generated in the 
cavity that in return may lead to entanglement between dif- 
ferent components of the system. In particular, stationary en- 
tanglement has been predicted between the cavity mode and 
a vibrating mirror [5-9], between an atomic ensemble or a 
Bose-Einstein condensate located inside an optical cavity and 
the vibrating mirror of the cavity [10-13], between two vibrat- 
ing mirrors of a ring cavity [14], between two dielectric mem- 
branes suspended inside a cavity [15], and between a mem- 
brane and a trapped atom both located inside a cavity [16- 
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18]. Further studies have addressed interesting problems of 
entangling mechanical oscillators [19], entangling optical and 
microwave cavity modes [20], and the creation of a photon by 
a vibrating mirror [21]. In this connection, we should men- 
tion the most recent work on the generation of entanglement 
in pulsed cavity optomechanics [22], and the work on the cre- 
ation of entanglement between two oscillating mirrors through 
the coupling of the mirrors to an atomic system [23]. 

The purpose of this paper is to explore coherence and cor- 
relation features of an optomechanical system composed of 
three bosonic modes realized with an one-dimensional opti- 
cal lattice located inside a single-mode cavity with a mov- 
able mirror. We are particularly interested in the effect of 
the first-order correlation (first-order interference) on entan- 
glement between two modes which is associated with second- 
order correlation functions. It is well known that in the para- 
metric down-conversion, served as a typical source of entan- 
glement, the signal and the idler beams are strongly entangled 
but always behave as mutually incoherent [24-26]. A similar 
conclusion applies to the correlations between modes of the 
optomechanical system, where it was demonstrated that the 
cavity and mechanical modes play the same role as the signal 
and the idler of a nondegenerate parametric oscillator and the 
modes behave as mutually incoherent [8, 11]. This seems to 
suggest that entanglement between two modes rules out the 
first-order coherence between them. 

We consider various situations where the modes of the op- 
tomechanical system can be made mutually coherent and to 
exhibit the first-order interference. This leads to an obvious 
question of to what extent of the first-order coherence could 
affect entanglement between the modes. To address this ques- 
tion, we use a polariton model of the optical lattice coupled to 
a single-mode cavity field and calculate various coherence and 
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correlation functions of the optomechanical system. We show 
that the system is capable of generating a wide class of coher- 
ence and correlation effects, ranging from the first-order co- 
herence, the anomalous autocorrelations and anomalous cross 
correlations between the modes. In a series of simple exam- 
ples we show that the generation of the first-order coherence 
between two modes of the system is equally effective in de- 
stroying entanglement between these modes. We illustrate our 
considerations by examining two and three mode interactions. 
After establishing the connection between the generation of 
entanglement and the first-order coherence, we consider the 
problem of the creation of entanglement between two inde- 
pendent modes by coupling them to an intermediate mode. 
We show that the coupling of the modes to the intermediate 
mode by a parametric interaction results in no entanglement 
between the modes. We then consider a different coupling 
configuration and find that an entanglement can be generated 
between two independent modes if one of the modes is cou- 
pled to the intermediate mode by the parametric interaction 
and the other is coupled by the linear-mixing interaction. 

The paper is organized as follows. We begin in Sec. II with 
the description of the model. We represent the finite size opti- 
cal lattice that is located inside a single-mode cavity in terms 
of Bloch-type waves called excitons, and diagonalize the in- 
teraction Hamiltonian of the excitons plus the cavity mode to 
describe the system in terms of bright polaritons. We then pro- 
ceed in Sec. Ill to study the dynamics of the system in terms 
of the quantum Langevin equations for relevant variables. We 
use the linearization approach to the equations of motion and 
arrive to a set of three couple differential equations for the 
fluctuation operators, which we solve for the steady-state. We 
apply the solution in Sec. IV to the calculation of the bipartite 
coherence and correlation functions of the polariton and the 
mechanical modes. We discuss the conditions for entangle- 
ment in terms of the squeezing fluctuations and the Cauchy- 
Schwartz inequality. The anomalous autocorrelation and cross 
correlation functions are introduced to discuss conditions for 
the violation of the Cauchy-Schwartz inequality. A possibility 
of generating a two-color entanglement is also discussed. In 
Sec. V we evaluate the first-order coherence and entanglement 
in the case of three coupled modes. Two coupling configura- 
tions of two independent modes to the intermediate mode are 
discussed. In Sec. VI, we examine parameter ranges in which 
the predicted coherence and correlation effects could be ob- 
served with the current experiments. Finally, in Sec. VII we 
summarize our results. 



II. THE MODEL 



We consider a finite size one-dimensional optical lattice 
located inside a single-mode cavity with one fixed partially 
transmitting mirror and one movable perfectly reflecting mir- 
ror, as shown in Fig. 1. The optical lattice is composed of TV 
regularly spaced and non-overlapping sites located at posi- 
tions r„ = nd, n — 1, . . . , TV, with total length L = Nd <C 
wo, where d is the separation between the sites and wo is the 



cavity mode waist at the position of the lattice [27]. The lat- 
tice is formed by two counterpropagating laser beams entering 
the cavity from the sides and forming a standing wave in the 
direction perpendicular to the direction of the cavity mode. 
The cold atoms loaded on the optical lattice are confined in 
an array of microscopic trapping potentials, forming a Mott- 
insulator-like medium with one atom per site [28, 29]. 

The atoms in the optical lattice are modeled as two-level 
systems with ground state \g n ) and excited state |e n ), sepa- 
rated by the transition frequency uj a and connected by a tran- 
sition dipole moment jl = (e n \f2 n \g n ), which can be assumed 
to be real with no loss of generality. Since the optical lat- 
tice is formed in the direction perpendicular to the cavity axis, 
the effect of the motion of the atoms (center of mass motion) 
on the coupling strength of the atoms to the cavity mode and 
on the radiation pressure on the movable mirror can be ig- 
nored. The situation would be different and the center of 
mass motion important if the optical lattice were generated 
along the cavity axis by the cavity mode [30], or by a stand- 
ing wave formed by running and reflected from the movable 
mirror laser beams [31]. 

The motion of the movable mirror is modeled as a quantum 
mechanical harmonic oscillator of mass m and resonant fre- 
quency ui m . The cavity mode is driven by an external laser 
field which is treated classically in our calculations and is 
characterized by its frequency lul and amplitude El- It has 
become common to consider the laser field as a source of the 
radiation pressure force on the movable mirror. 




FIG. 1. (Color online) Schematic diagram of the system. An optical 
lattice composed of regularly spaced atoms is located inside a single- 
mode cavity driven by a laser field. The cavity is composed of one 
fixed and one movable mirror that can undergo harmonic oscillations 
due to the radiation pressure induced by the laser field. The atomic 
sites have parallel dipole moments /2 oriented in the direction making 
the angle a with the lattice direction. 

The total Hamiltonian H of the system can be written as 

H = H c + H ex +H + H I , (1) 

where 

H c = Hu) c a'a (2) 
is the free Hamiltonian of the cavity mode 

H ex =hJ2 u a BlB n + J a B\B m , (3) 
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is the Hamiltonian of the electronic excitation in the atoms of 
the optical lattice 



where 



H 



1 



huj m (q 2 + p 2 ) , 



(4) 



is the free Hamiltonian of the oscillating mirror, and 

Hj = h'y^gn (B n a + a* ~ hGoo'aq 



h{E L a)c~ luJLt + E* L ac luJLt ) 



(5) 



is the interaction Hamiltonian of the cavity mode with the 
electronic excitations, the movable mirror and the external 
laser field, respectively. 

Here, a^ and a are creation and annihilation operators of 
the cavity mode of frequency lj c . The operators q and p are, 
respectively, the dimensionless position and momentum op- 
erators of the oscillating mirror that satisfy the fundamental 
commutation relation [q,p] = i. The contribution of elec- 
tronic excitations in atoms is expressed as a sum over normal 
Boson creation and annihilation operators, B n and B n , re- 
spectively, one for each site n of energy tiui a , and 
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, 1 — 3 cos a) 



(6) 



represents the contribution of the nearest-neighbor dipole- 
dipole interaction between atomic sites, with parallel dipole 
moments oriented in the direction making the angle a with 
the lattice direction. In a practical situation, the angle a is 
fixed by the polarization direction of the cavity mode. At low 
number of electronic excitations, that we consider here, the 
operators B n can be treated as bosonic operators. 

The first term in the Hamiltonian (5) describes the inter- 
action of the electronic excitations with the cavity field. The 
strength of the interaction is characterized by the Rabi fre- 
quency g n . The interaction retains only the terms which play 
a dominant role in the electric-dipole and rotating-wave ap- 
proximations. The higher order and antiresonant terms which 
would make much smaller contributions have been omitted. 

The second term in the Hamiltonian (5) describes the op- 
tomechanical radiation-pressure interaction which couples the 
cavity photon number n c = a) a to the position operator q 
of the oscillating mirror with the coupling strength Go = 
(uj c / L c ) \Jh,j muj m , where m is the effective mass of the me- 
chanical mode and L c is the length of the cavity. Finally, the 
parameter El describes the coupling strength of the laser field 
to the cavity mode. 

We now consider the energy states of the optical lattice that 
is determined by the Hamiltonian (3). It is easily verified that 
in the bare basis of the lattice sites, {|£/i), ■ • ■ , \Un)}, 
where \Ui) = \a) Wj^i \9j), the Hamiltonian (3) is not diag- 
onal due to the presence of the dipole-dipole interaction J a . 
The diagonalization of the Hamiltonian (3) results in eigen- 
states described by Bloch-type waves, called the collective ex- 
citation modes or shortly excitons. The diagonal Hamiltonian 
is of the form 



Hex — ^fcCfcCfc, 



(7) 



u k =0J a + 2J a cos 



■nk 
N+ 1 



(8) 



is the frequency of the fcth exciton mode, and C k and Ck are, 
respectively the creation and annihilation operators of the ex- 
citons. The operators C\ and Ck are obtained from the cre- 
ation and annihilation operators of an electronic excitation at 
site n using the transformation 



B n 



N 



1 ^ 

k 



imk 
N+ 1 



C k . 



(9) 



Simply, the Ck operators are obtained by inverting the above 
transformation. 

Thus, in terms of the exciton operators, the total Hamilto- 
nian of the system can be written in the form 



H =i? + if/, 
where Ho is the free Hamiltonian 

Ho = hui c a^a + hu k C\Ck + ^huj m (q 2 



(10) 



(11) 



and Hi is the interaction Hamiltonian 
H I = hY J fk (C k a) + H.c. 



hGoa) aq 



odd k 



with 



fk 



h (E L a)e- %UJLt 



E* L ae iuiLt ) 



he V(N + 1) 



■ cot 



71 k 



2(N+1) 



(12) 



(13) 



Note that the interaction involves exciton modes with odd k 
only [27]. In addition, the coupling constants fk are not iden- 
tical, so that the exciton modes k are not equally coupled to the 
cavity mode. Since the cotangent function decreases rapidly 
with k, one can easily verified that the strength of the coupling 
of the k = 1 mode is stronger by k 2 from the k ^ 1 modes. 
This indicates that only the k = 1 mode can be strongly cou- 
pled to the cavity field, with the modes k ^ 1 weakly coupled 
to the field. 

The strong coupling of the k = 1 exciton to the cavity mode 
prompts us to write the Hamiltonian in the form 



H — Hi + H2 , 

where 

H x = Ku c a*a + hw^ld + hfi (CV + H.c.) 

and 



(14) 



(15) 



H-2 = —huj 



p 



2 ) - hGoah 



oa' aq 
z, 

h{E L a^c-^ Lt + E* L ac tuJLt ) . 



(16) 
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We may diagonalize the Hamiltonian Hi to find new operators 
of the combined k = 1 exciton plus the cavity field system. 
The cavity mode can be considered to "dress" the exciton and 
to form along with it a single "polariton" quantum system. 
This reflects the fact that photons are exchanged between the 
exciton and cavity modes. The dressed operators, the eigen- 
operators of the Hamiltonian Hi , are found by the following 
unitary transformation 

4< = (cos0)Ci — (sin cj))a, 

$ = (sin0)Ci + (cos (j))a, (17) 

where the rotation angle <f> is defined by 

cos^=i + A (18) 

with 5 = (cj c - Wi)/2 and fl = (jf + <S 2 )i The angle cj) 
belongs to the interval [0,7r/2]. The polaritons are coherent 
superpositions of the exciton and the cavity fields. For 5 = 0, 
these are equally weighted, maximally entangled, superpo- 
sitions of the fields, whereas for a large positive detuning, 
6 ^ /i> the exciton and the cavity field disentangle that then 
then polariton \I/ becomes purely atomic (excitonic), while the 
polariton $ becomes purely photonic. 

In terms of the polariton creation and annihilation opera- 
tors, the Hamiltonian ( 14) takes the form 

h = h(uj - + h(uj + 

+ - f^ m (q 2 +P 2 ) - hGoalaq 

+ h (^aV^ 1 * + E* L ac tulLt ) , (19) 

where ujq = (uj c + lo\)/2 is the mid-frequency of the two 
polariton modes, and the annihilation operator for the cavity 
mode is related to the annihilation operators of the polariton 
modes by a = $ cos <fi — sin <j>. 

This shows the familiar coupled exciton-photon mode split- 
ting [33]. The mid-frequency of the polariton modes is an 
average value of the cavity frequency and the excitonic fre- 
quency. Thus, wo is always between tu c and tui, that the fre- 
quency ojq is pulled away from the cavity frequency towards 
the excitonic frequency. When the cavity frequency is tuned 
to exact resonance with the excitonic frequency, i.e., u> c = u>i, 
there is no mode pulling, i.e. ljq = oj c . 

III. LINEARIZED FLUCTUATION ANALYSIS 



Given the Hamiltonian of the system, we now proceed to 
study the dynamics of the system in terms of the Heisen- 
berg equations of motion for relevant variables, the polari- 
ton modes and the mirror mode operators. However, a proper 
analysis of the system must include losses due to the coupling 
of the cavity mode, the polaritons and the oscillating mirror 
to their local environments. Therefore, we introduce so-called 
phenomenological damping terms that the cavity field ampli- 
tude is damped with a rate k, the amplitude of the exciton's 



field is damped with the atomic spontaneous emission rate j a , 
and the oscillations of the cavity mirror are affected by the 
quantum Brownian noise acting on the mirror leading to the 
damping of its oscillations with a rate j m . The inclusion of 
the losses to the Heisenberg equations of motion results in a 
set of nonlinear Langevin equations, that written in the frame 
rotating at the laser frequency ujl have the form 

q = U m p, 

P = -7mP - u m q + G s ^y + G c $ f $ 

* = -E* L sin0 - [7 + « (A L - 0)] * + iG s g* 

<j> = E* L cos - [7 + i (A L + 0)] $ + iG c q$ 

- ^iGqV + v^$ m , (20) 

along with the corresponding equations for the adjoint opera- 
tors & and Here, A L = uj -uj l , G s = G sin 2 0, G c = 
Go cos 2 4>, G = Go sin (20), and we have chosen j a = k = 7. 
Moreover, the laser field is assumed to be tuned close to 
the cavity and atomic resonance, in the sense that the detun- 
ing Al is small compared to the optical frequency us^. We 
shall assume additionally that the input noises to the polariton 
modes, and which are the sum of the input noises to 
the cavity and to the excition modes, are frequency indepen- 
dent Gaussian (white) vacuum noises, so that all first moments 
vanish, (* m ) = (tf J n ) = ($ m ) = = 0, and only non- 

zero are the following second moments 

(*in(ml(t')) = (*in(ml(t')) = S(t - t 1 ). (21) 

Similarly, for the quantum Brownian noise £(i) which arises 
from the coupling of the oscillating mirror to its local envi- 
ronment, we assume that it is a frequency independent Gaus- 
sian thermal white noise, so that all first moments vanish, 
= 0, and the only non-zero are the following second 
moments 

\ mmn + (mm)) = ( n +l) w - o. (22) 

where n = [exp(fiw m /fcsT) — is the mean number 
of the thermal excitations at the frequency of the mechanical 
mode, ks is the Boltzmann constant and T is the temperature 
of the environment. 

The polaritons ^ and $ might reasonably be called "bright" 
polaritons since they are damped with the rate 7. It is clear 
from Eq. (20) that in the absence of the mechanical oscillator 
the system would consist of two completely decoupled bright 
polaritons. The effect of the mechanical oscillator that inter- 
ests us most here is to introduce both shifts of the frequencies 
and coupling between the polaritons. 

The exact treatment of the problem that involves quantum 
properties of the mechanical oscillator requires to deal with 
the system of nonlinear differential equations. The system 
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of the equations is difficult to solve. Therefore, we use the 
linearization approach [32] by assuming that each operator of 
the system can be written as the sum of its steady-state mean 
value and a small fluctuation around the steady-state 



q = q s + Sq, p = p s + 5p, 



(23) 



Note that this linearization approach is equivalent to the as- 
sumption of Gaussian distributions that describe fluctuations 
of the system around its stationary state. In this approach, 
Eq. (20) decouples into a set of non-linear equations for the 
steady-state values and a set of differential equations for the 
fluctuation operators. 

We first determine the average steady-state values of the 
operators. By taking the mean values of the operators, and 
then by setting the left-hand sides of Eq. (20) to zero, we ob- 
tain the following steady-state solutions for the oscillator vari- 
ables 

- ^9 



Ps 



0. 



ty s sin < 



(24) 



and the steady-state values of the polariton fields are found 
from the solution of two coupled nonlinear equations 



E* L sin 4> = 


7 + 


i (a, 


- n\ 


* s 


E* L cos 4> = 


7 + 


*(a, 


+ n) 





iGq s $> s , 
iGq s V s , (25) 



where A q = Al — ^q s Go and Q = fi — ^q s Go cos(20). 

Under the linearization procedure and introduce the annihi- 
lation operator of the oscillating mirror, Sb= (5q + i5p)/V2, 
the Langevin equations for the fluctuation operators satisfy the 
following set of differential equations 

Sb = - Qt™ + iojr^j Sb + ^ m 5tf - ^iG* (<f¥ + 5^) 

+ ^'G$ (<5$ + S&) + Vt^£, 
6$ = - (7 + »A*) 6® - ^G* (5b + Stf) 



iG q S<I> 



<5$ = 



iG q 6V 



2t~$ 4 



- (7 + iA») <5$ + -%G 9 (6b + 6b*) 



(26) 



where we have chosen A* = (A q — Cl) and A$ = 
(A q + fi), and have introduced the abbreviations G<j = 
n/2(|G$ s -G s * s ), G* = x/2(G c $ s -iG* s ) and G q = 
Gq s /2. Then we can get G$ = G$tan</>. Also we see 
that G* and G$ are the effective coupling constants of the 
polaritons \I> and $ to the cavity field, respectively. 



IV. NANO-MECHANICAL ENTANGLEMENT 



We now apply Eqs. (26) explicitly to search for entangle- 
ment and correlations between the modes. Notice the presence 
of three different frequencies at which the fluctuation opera- 
tors oscillate, u> m , Ay and A$. Thus, depending on whether 
we would like to entangle one or both polaritons to the oscil- 
lating mirror, we should choose w m to match to either the fre- 
quency of one of the polaritons or to the frequency A q , which 
is the mid-frequency of the two polaritons. To illustrate this, 
we introduce a rotating frame through the relations 

Sb = Sbe lUm \ = <5#e iA **, 5$ = 6^c iA ^, (27) 

and find that in the rotating frame Eqs. (26) become 

Sb = - ^ lm 6b + ^ lm 5bh 2 ^ 1 + VT^e^' 

- -iG*^#e-^ A *- Wm) * + s¥c l ^ +UJ ^ 

+ Ug^^C- 1 ^-^' + J$t e i(A*+a, m )tj ^ 

5* = - 7<5§ - iiG* (Sbe^*-"™^ + <ft,V(A*+« m )t) 

- iG 9 (5$e 4(A *" A * ) * + v^y^^e^**, 

5<t = - 7 <5$ + -iG$ (sbe i{A *- Um)t + stf c i{A * +UJm)t ) 



iG (? J*e l(A *^ A * )t + V27*me l 



(28) 



We see that the coupling terms of the polaritons 8^f and <5$ to 
the mirror operator oscillate in time with frequencies A$ ±w m 
and A$ ± us m , respectively. When the equations are inte- 
grated over any measurable time interval, these oscillatory 
terms make a negligible contribution and can be ignored if 
they are different from zero. Therefore, in order for the cou- 
pling effects to be significant, we must have A^ = ±oj m or 
A$ = ±Lu m , when either of these terms become independent 
of time. These are optimal conditions for coupling of the mir- 
ror to either iff or $ polariton. 

An alternative choice of the new rotating frame 



6b = Sbe iUmt , 8¥ = 5* 8$ = <5$ e iA "*, 
results in the following transformed equations 

Sb = - ^ m 8b + ^mSbU 2 ^ + Vr^e 1 ^* 



(29) 



2 
1 

+ 2 
<5# = - ( 7 - 



5<f> 



(7 - iQj 5$ - iG q 5$ + y/¥f^ in e iA ^ 
^iG* (Sbe 1 ^-^ + 8PS A " +Um)t 
(7 + in) 6$ - iG q 5V + v/2^$ m e' 
izG$ (<56e i ( A »- w '")* + 5St e »-(A,-H-m)t) 



JA„t 



(30) 



Now the exponential factors have frequency centered on A q , 
and in marked contrast to the previous situation, the choice of 
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either A q = ui„ 



or A q = —ui m would result in the simulta- 



neous coupling of both polaritons to the oscillating mirror. It 
is interesting to note that the coupling between the polaritons 
is independent of the choice of the frequency u) m . 



A. Bipartite polariton-mirror coupling 

Let us first examine a bipartite coupling between the polari- 
ton "J and the oscillating mirror. According to Eq. (28), this 
could be achieved by choosing A* = — u) m , which has been 
shown as a necessary condition for entanglement between two 
bosonic modes [8]. In this case, the two modes are coupled 
by a parametric interaction [34, 35]. The tuning of A* = io m 
would not produce entanglement between the mirror and the 
polariton modes, since in this case the two modes are coupled 
by a linear-mixing interaction. 

By choosing A^ = — u> m , after discarding the rapidly os- 
cillating terms, we find from Eqs. (28) that the set of the 
differential equations (28) can be simplified to two separate 
sets of coupled differential equations for pairs (5b' , 8^f) and 
(8b, 8^'). The equations of motion for the pair (5w , 8^f) are 



<5* 



27?(«), 

27*in(t), (3D 



where = ^exp(-iu m t), = * m cxp(-iw m i), 

and we have put j m = 27. Note that the two coupled modes 
oscillate at the same frequencies. This may result in the so- 
called one-colour entanglement, i.e. entanglement between 
two modes of the same frequency. 



B. Identification of entanglement from the squeezing condition 

In order to show that entanglement and squeezing can be 
created between the ^ and b modes, we introduce in-phase 
and out-of-phase quadrature components of the fluctuation 
operators 



5V X 



<5* 



V V2i 
1 



8$ x = -^(S^e-^ + S^e** 
72V 

8$ y = JL^e-** - $*t e *V^ 



(32) 



and 



8A X = ±= (s$ x - 5q) , 8T X = -L (s¥ x + Sq 



(33) 



In order to see if there are correlations existing between 
the polariton and the mechanical mode b, we must examine 
properties of the sum operators, either SY X or 8T y and the 
difference operators, either 5A X or 8A y that act on both sys- 
tems. The variances in the sum and difference operators are 
given by 



A(8T 
1 

~~ 2 



= (STl) - (8T X 

A(8V x ) 2 + A(Sq) 2 



2((8^ x 8q)-(8^ x )(8q) 



A (6 Ax) = (SAl) — (SA X ) 2 

= I \a(6V x ) 2 + A(Sq) 2 - 2 ( (5V x 8q) - (SV x )(8q) 



(34) 



The product of these two variances is 

A^T^'A^A,) 2 -^ [A(^) 2 



(6* x 8q) 



+ A(Sq) 2 



Hence, we have an uncertainty relation 

A(S* x f + A(Sq) 2 



A(Sr x ) 2 A(SA x ) 2 < x - r 



(35) 



, (36) 



with equality only holding if the systems are uncorrelated. 
Solving Eq. (31) for the steady-state (t — > 00), we find 



and 



A(ST x y 



A(SA X 



7(n 


+ 1) 


2(7- 


\G^) 


7(n 


+ 1) 


2(7 4 




7 2 (n 


+ 1) 


( 7 2_ 




7 2 (n 


+ 1) 



(37) 



(7 2 - \G%) 



(38) 



It is seen that the variances A(8^> x ) 2 and A (8q) 2 are both 
larger than their vacuum limits, i.e. 1/2. Thus, both modes 
are not themselves squeezed and display thermal fluctuations. 
However, substituting Eqs. (37) and (38) into Eq. (36), we 
find that the polariton and the mechanical modes are corre- 
lated when 



7 



> 1, 



(39) 



which is always satisfied as long as 7^ 0. We stress that 
the above inequality is necessary but not sufficient condition 
for squeezing (entanglement) between the modes. In other 
words, the modes could be correlated but not enough to beat 
the quantum limit for fluctuations. Equivalently, we may say 
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that the modes are correlated classically and the quantum limit 
can be beaten only if the modes exhibit quantum correlations. 

The sufficient condition for squeezing is that the variance 
A (5T X ) 2 is reduced below the limit for quantum fluctuations, 
i.e. below 1/2. It is easily verified from Eq. (37) that the cor- 
relations between the modes will lead to squeezing in the su- 
perposition ST X when > 2n-f. However, there is an upper 
limit on Gi> imposed by the condition of stable steady-state 
solutions ofEqs. (31). One can easily find from Eqs. (31) that 
the stability of the steady-state solutions requires G^ < 2-f. 
Thus, combining the stability and squeezing conditions, we 
find that the necessary and sufficient condition for entangle- 
ment between the polariton ^ and the vibrating mirror mode is 



2n-f < G* < 2j, 



(40) 



which, on the other hand, indicates that the modes can be en- 
tangled only if n < 1. 

Since the modes are correlated for any < 2j and the 
condition for squeezing (entanglement) is that G* must be 
greater than 277,7, there is evidently a significant restriction 
on the strength of the coupling of the vibrating mirror to the 
polariton mode. The condition (40) for entanglement is es- 
sentially similar to that of a microcavity mode and a vibrating 
mirror treated by Vitali et al. [8]. 



C. Violation of the Cauchy-Schwartz inequality and 
anomalous correlations 

An alternative and in fact more elegant way to study entan- 
glement between two modes (A, B) is the Cauchy-Schwartz 
inequality [36] 



X(A,B) 



(2) (2) 

9 A 9b 



(2) 

9ab 



> 1. 



(41) 



Here, X(A,B) is the so-called Cauchy-Schwartz parameter. 

(A^B^AB) 



(2) 
9ab 



(42) 



(A^A)(B^B) 

is the normalized second-order cross correlation function, and 

(A^ 2 A 2 ) (2) (BVB*) 
9b - 



(2) 
9 A 



(A1A)' 



(B*B) S 



(43) 



are the normalized intensity autocorrelation functions of the 
modes A and B, respectively. 

Since the modes obey the Gaussian statistics, the correla- 
tion functions can be readily related to coherence functions 



9 (2) 



\V(n,n)\ 



9AB = 1 + h(A,B)\ 2 + \V(A,B) 



A,B, 
2 



where 7(a,b) > defined as 



|(AtB)| 



(44) 



(45) 



is the degree of the first-order coherence, 

, \(A 2 )\ , , \(B 2 ) 



\V(A,A)\ 



(A^ A)' |7/(S ' B)I (BIB)' 
are degrees of the so-called "anomalous" autocorrelation, and 



(46) 



\V(A,B) 



\(AB)\ 



(47) 



is the degree of the anomalous cross correlation [37-42]. 
Equation (44) shows that the second-order autocorrelation 
functions depend on the anomalous autocorrelation whereas 
the second-order cross correlation function depends on the 
first-order coherence and the anomalous cross correlation. 

We now determine conditions under which the correlation 
between the polariton and the mechanical mode b violate 
the Cauchy-Schwartz inequality (41). Following the method 
introduced in the book [43], we can easily solve Eqs. (31) to 
find that in the steady-state 



(<5* 2 ) = (5b 2 ) = 0, (5^ 5b) = 0, 



and 



(5^ 5b) = --(n 



1) 



7G* 



4 U * 



(48) 



(49) 



which together with 



(5&5V) = - (n + 1) 
8 

(SbUb) = \ 



(n-l) + 



ir*2 ' 
(n + l)7 2 



r 



1^2 



(50) 



give the following condition for the violation of the Cauchy- 
Schwartz inequality 



(?i + l)7 2 



2n7 2 



> 1. 



(51) 



We emphasize that the inequality (51) is the necessary and 
sufficient condition for entanglement between the modes. It 
is easily verified that the inequality is satisfied as long as 
n < 1 and G* < 27. This shows that the Cauchy-Schwartz 
inequality is violated under the same conditions the squeez- 
ing is generated between the modes, see Eq. (40). Moreover, 
the inequality (51) implies that the violation of the Cauchy- 
Schwartz inequality is achieved simply by the requirement 
that the anomalous cross coherence is greater than the product 
of the intensities of the modes. It is interesting to note that the 

condition (51) corresponds to the case of each of the modes 

(2) (2) 

being in the thermal state, g^, ' = <?> ' = 2, with no the first- 
order coherence between them, 7(^^)1 = 0. Thus, the fields 
of the polariton and mechanical modes are entangled but there 
is no interference. 

We now illustrate the above predicted limitations for the 
occurrence of entanglement and also provide the quantita- 
tive value of the bipartite entanglement. For this, we use 
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the logarithmic negativity that is known as the necessary and 
sufficient condition for entanglement of two-mode Gaussian 
states [44, 45] 



E N =max{0,-log 2 [2V S ]} , 



(52) 



where V s is the smallest sympletic eigenvalue of the partially 
transposed correlation (covariance) matrix V, with elements 



Vij = (ui(oo)itj(oo) + Uj(oo)iij(oo))/2, 



(53) 



where u;(oo) is the steady-state value of the ith component of 
the vector u: 



8q,Sp, Sy x ,S^, 



(54) 



with Sq = (5b + S~tf)/V2,5p = (5b - 5#)/V2i, and 
5^> x , S^/y defined in Eq. (32). 

The steady-state values are readily calculated using the 
equations of motion (31), from which we get the following 
matrix equation 



(55) 



u(t) = Au(t) + V2 7 7f(i), 








is given 


by 










— ^Gty 


-7 - 









-7 











-7 



(56) 



and 



m = Uin(t),Pin(t),^Ut)^ y in (t) 



(57) 



The matrix equation (55) is a simple first order differential 
equation with time-independent coefficients, and is solved by 
a direct integration. The formal solution is given by 



u(t) = u(0)e At + V27 / dt' f)(t - t') 



(58) 



where u(0) is the vector of initial values of the components. 
For the steady-state, we take the limit of Eq. (58) as t — >• oo. 

Since the noise £(t) is (^-correlated, so that it describes a 
Markovian process, the steady-state correlation matrix is then 
derived from the following equation [32]: 



AV + VA J 



D. 



(59) 



where D = diag[(2n + 1)7, (2n + 1)7, 7, 7] is the diffusion 
matrix stemming from the noise correlations. 

Figure 2 illustrates the dependence of the logarithmic nega- 
tivity on n and G*. It is apparent that the modes are entangled 
for n < 1 and < 27. Once n is greater than 1, entangle- 
ment becomes impossible irrespective of G^. 

The dependence of the second-order correlation func- 
tions (44) on different kind of coherence functions allows us 
to determine which of the coherences work for and which 
work against the creation of a strong entanglement between 




G h 



FIG. 2. 

and G* /j. 



Variation of the logarithmic negativity Em with n 



the modes. A simple analysis of Eqs. (41) and (44) shows that 
in general the Cauchy-Schwartz inequality is violated when 

(2+\V(A,A)\ 2 )C2+hB,B)\ 2 ) < (1 + |7(A,B)| 2 + I^,B)| 2 ) 2 - 

(60) 

It follows from Eq. (60) that the general condition for the vio- 
lation of the Cauchy-Schwartz inequality cannot be viewed as 
exclusively dependent on \t](a,b) | 2 - One could notice that the 
Cauchy-Schwartz inequality would be easier to violate if the 
first and second order cross correlations were simultaneously 
created while the anomalous autocorrelations were kept zero. 
However, we shall demonstrate below that this situation is un- 
founded, that the simultaneous creation of the cross correla- 
tion functions is equally effective in creating the anomalous 
autocorrelations. 



D. Squeezed vacuum environment 

When evaluating the correlation function (6wSb), 
one can find that the function depends solely on the 
noise two-photon correlation functions (^l n (t)^l n (t')) and 
(I* (*)£*(*')■ Therefore, the correlation function (S^Sb) 
could be different from zero when either the polariton or the 
mirror were located in the environment whose modes exhibit 
nonzero two-photon correlations. An example of such envi- 
ronment is a squeezed vacuum field which in the case of the 
mirror environment is determined by Eq. (22) and the follow- 
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ing second moments 

(i(t)iHt')) = (n+l) S(t - t'), <f (*)£(*')> = nS(t 0. 
(£(t)£k') = mS(t + if), ( t )£t(f ) = m *A(f + f), (6i) 

where n is the squeezing photon number and |m| < 
\Jn(n + 1) measures the strength of two-photon correla- 
tions [46]. Thus, if the mirror were oscillating in the squeezed 
vacuum field, this would create the first-order coherence be- 
tween the modes. 

Let us apply these considerations explicitly to the first- 
order correlation function (S^Sb). If the mirror oscillates 
in the squeezed vacuum field, we readily find that 



(6W6b) 



t AM - 



(62) 



which is nonzero as long as the polariton and the mirror are 
coupled to each other. When Eq. (62) is substituted into 
Eq. (45), we find for the degree of the first-order coherence 



|7(*,i>)| 



7|m| 



{(n+1) [2n^-\{n-l)G%]} 



1/2 ' 



(63) 



which is less than 1 in general, and becomes unity only in 
the limit of G* = 2j and fi > 1, Expression (63) shows 
that the interaction of the oscillating mirror with the squeezed 
vacuum field results in the first-order coherence between the 
polariton and the mechanical modes. That happens because 
the two-photon correlations present in the squeezed field have 
the effect of inducing stimulated two-photon processes which, 
it turns out, are sufficient for the polariton and mirror fields to 
become mutually coherent. Equivalently, the oscillating mir- 
ror that scatters photons from the squeezed vacuum to the po- 
lariton mode gives rise to phase locking between the polari- 
ton and mechanical modes. The correlations in the squeezed 
vacuum field have therefore induced coherence between the 
polariton and mechanical modes. It is interesting to note that 
the squeezed correlations do not effect the anomalous cross 
correlation (S^Sb). 

The creation of the first-order coherence should, according 
to Eq. (60), enhance the violation of the Cauchy-Schwartz in- 
equality. However, this is not the case, because the squeezed 
vacuum not only creates the first-order coherence but also the 
anomalous autocorrelations. It is easy to find that 



(^ 2 > = 



(5b 2 ) = -m 



(7 2 - \G%) 



(64) 



which shows that in both modes the anomalous autocorrela- 
tions are induced by the squeezed vacuum. The degrees of the 
anomalous coherences are then given by 



I^C*,*) I 
\V(b,b)\ 



m 2 )\ = b 

(<J#t<S*) n 
\(8b 2 )\ _ |i 



(27 2 



jG%) 



(SbUb) [2n 7 2 -±(n-l)G|] 



(65) 



Clearly the squeezed vacuum conspires to create the anoma- 
lous autocorrelation functions. Note that |7(*.&)| < 

\/|??(*.*)||?/(b,6)| and the equality holds at the threshold value 
of G^, = 27. It can also be seen that the anomalous auto- 
correlation of the mechanical mode is greater than that of the 
polariton mode and the equality between the autocorrelations 
is achieved when G* = 27. 

The presence of the anomalous autocorrelations alters the 
condition (60) for the violation of the Cauchy-Schwartz in- 
equality. It is clear that the left-hand-side of Eq. (60) is en- 
hanced by the anomalous autocorrelation functions. Thus the 
inequality (60) can be harder to achieve. Needless to say, 
the creation of the first-order coherence among the modes is 
achieved at the cost of a corresponding decrease in the entan- 
glement between the modes. 




FIG. 3. Variation of the Cauchy-Schwartz parameter X(*,6) w i m n 
and G* / 7 for \m\ = y/n(n + 1). 

Figure 3 illustrates the Cauchy-Schwartz parameter as a 
function of n and Gj, for the maximally correlated squeezed 
vacuum field, \m\ = \Jn(n + 1). The effect of including 
the squeezing correlations is clearly to restricts the range of n 
and G* over which the modes are entangled. 

We may conclude this section that for the best conditions 
to entangle two degenerate modes is a situation of mutually 
incoherent modes each being in the thermal state. 



V. THREE-MODE COHERENCE AND ENTANGLEMENT 



In this section we determine conditions for correlations and 
entanglement when both polaritons ^ and $ are simultane- 
ously coupled to the mechanical mode. This is a situation of 
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the three-mode interaction and appears in the system when the 
frequency of the mechanical mode w m is tuned to the midd- 
frequency of the two polaritons, i.e. when uj m = — A q . We 
also briefly study the two-colour entanglement that may occur 
between two modes of different frequencies. 

When we chose A 9 = — ui m and make the secular approx- 
imation, we find that Eqs. (30) simplify to two separate sets 
of three coupled differential equations for (Sw , 5^, 5$) and 
(5b, 8^', 5<&'). For example, the equations of motion for the 
set (5b\5^,5^) are 

6& = -p m S# + \iG^,S^> - ^G $ <5$ + VT^fW, 

6$ = - (7 - i&) ^ - iGqS® - -iG 9 6# + ^y^ in (t), 

5$ = - (7 + itlj 5$ - iG q 6$> + -iG&Stf + ^y^ in (t), 

(66) 

where 3>j„(t) = <E>m exp(— iu> m t). Equations (66) are quite 
different from Eq. (31) that now both polaritons are coupled 
to the mechanical mode and the coupling is of the type of a 
parametric interaction. Normally, we would expect that this 
kind of coupling should result in an entanglement between 
the polaritons. As we shall see below, this hope is unfounded, 
the coupling of the polaritons with the mechanical mode by 
the parametric interaction results in the first-order rather than 
a second-order coherence between the polaritons. 

Before going into detailed studies of the conditions for en- 
tanglement between the modes, we first comment about cer- 
tain general features of the simultaneous coupling of the po- 
laritons to the mechanical mode that follow from Eq. (66). We 
see that the effect of the mechanical mode on the dynamics of 
the polaritons is twofold. Firstly, the mechanical mode cou- 
ples the polaritons to each other with the coupling strength G q . 
This indicates that in the presence of the oscillating mirror, 
the polaritons \I/ and $ are no longer the eigenstates of the 
system. Secondly, the mirror couples to the polaritons with 
different coupling strengths G* and G$. The magnitude of 
the coupling strengths depends on whether a given polariton 
is maximally entangled or not. When the polaritons are max- 
imally entangled G* = G$, otherwise G<j ^ G$ when the 
polaritons are non-maximally entangled. 

The presence of the coupling between the polaritons 
prompts us to make an unitary transformation to obtain 'new' 
orthogonal polariton modes. It is easily shown that the anni- 
hilation operators of the orthogonal superposition modes are 
of the form 

50 = cos(0 + tp)8C\ — sin(0 + <p)5a, 

511 = sm((j> + (p)Sd + cos(> + ip)Sa. (67) 

where the angle ip is defined by 

, 1 n 

cos^ = - + -, (68) 

with U = \J& 2 + G^. Note that the angle ip belongs to the 
interval [0, 7r/4]. 



Several interesting features can be found in Eq. (67). 
Firstly, we note that the initially maximally entangled exci- 
ton and the cavity fields, i.e. <j> = it/ A and G q = 0, become 
non-maximally entangled when the mechanical mode is in- 
cluded, G q 7^ 0. Thus, the effect of the mechanical mode on 
the initial maximally entangled exciton and the cavity modes 
is to destroy the superposition (polariton) modes. Secondly, 
if initially the exciton and the cavity fields were disentangled, 
i.e. = 0, they remain disentangled even in the presence of 
the mechanical mode because in this case G q = and then 
ip = 0. Thirdly, an initial non-maximally entangled state be- 
tween the exciton and the cavity modes, < <j> < ir/2 can be 
transferred to the maximally entangled state by the mechan- 
ical effect. It happens when <fi + p = it /A. However, the 
maximally entangled state can be reached only for initial su- 
perpositions with <j> < 7r/4. Otherwise, for <f> > it /4, the 
superposition cannot be transferred by the mechanical effect 
into the maximally entangled state. 

We now turn to the calculation of the correlation functions 
between the modes which provide the information about co- 
herence and entanglement between the modes. In terms of the 
transformed operators 50 and SH, the equations of motion 
Eqs. (66) take the form 

( 5tf\ ( 1 -\iG e -±iG v \ /5tf\ 

5Q = - hG e fr + iU) 56 
\5tlJ V ¥ G - (1-iU) J \5llJ 

+ V^y O m (t) , (69) 

\n m (t) / 

where Gg = (G* cos ip + Gj> sin ip) and G^ = (G* sin <p — 
G$ cos tp) are effective couplings of the transformed polari- 
ton modes to the cavity mode. It is clear from Eq. (69) that 
the transformed polaritons are independent of each other, they 
oscillate with different frequencies, shifted from u> m by an 
amount ±J7, and both are coupled to the mechanical mode 
by the parametric amplification process. Notice that Gg = 
— Gtt tan(0 + ip). 



A. Two-colour entanglement 

Since the modes oscillate at different frequencies, we may 
apply the equations (69) to the problem of two-colour entan- 
glement [47-53]. If we assume that G^ = 0, we then can ig- 
nore coupling of the 5H polariton to the mechanical mode and 
limit ourselves to considering the case of two coupled modes 
of different frequencies. This could takes place, for exam- 
ple, when the cavity field frequency were on resonance with 
the exciton frequency ((5 = 0) and the coupling of the mir- 
ror to the cavity mode were much stronger than the coupling 
of the atoms to the cavity mode, i.e. G q 3> In this case, 
<f) = ip = it/4 and then it is easy to find that G n = 0. 

To check the conditions for the stable steady-state solu- 
tions, we put G-k = in the 3x3 matrix appearing in Eq. (69), 



11 



and find the following eigenvalues 



A x = 7 - iU, A 2 , 3 = (7 ± V G e - u * J + ( 7 °) 

We see that a threshold occurs for the coupling strength at 
Gg = U where the eigenvalues A2,3 change character. Be- 
low the threshold, Gg < U, the real parts of the eigenval- 
ues are positive irrespective of the values of the parameters 
involved. Above the threshold, Gg > U, the real parts are 
positive when \jG 2 e — U 2 < 7. We should stress here that 
under the condition of G q ^ Cl, taken here, the inequality 
Gg < U, i.e. the below threshold situation always holds. 
Therefore, as long as Gg < U there are no restrictions on the 
parameters for the the system to decay into a stable station- 
ary state. However, we have seen in Sec. IV that in the case 
of two degenerate modes there was a severe restriction that a 
stable steady state only exists for very weak optomechanical 
couplings, Gg = \f2Gq, < 2^/2j. 

After establishing the stability conditions, we now solve 
the set of the resulting equations of motion for the steady-state 
and find the following non-zero correlation functions 



(StfSb) = n 
(SbSO) = - 



(n + l)Gg 



2(U 2 -G 2 + 4 7 2 )' 

1 (n + 1) (27 - iU)G 6 

2 (C/ 2 -G 2 + 4 7 2 ) : 
(n + l)G 2 e 



2 (U 2 - G 2 + 4 7 2 



(71) 



If we now substitute Eq. (71) into Eq. (47), we arrive at 
the following expression for the anomalous cross correlation 
function 



\V(b,e) 



{n + 1) (4 7 2 + U 2 ) 
[2n (4 7 2 + U 2 ) -{ft -I) Gg 



(72) 



The expression (72) differs markedly from the one we encoun- 
tered in Eq. (51), although they become identical in the limit 
of degenerate modes, when U = 0. Similar to the case of de- 
generate modes, the condition for the violation of the Cauchy- 
Schwartz inequality, |?7(b.e)| 2 > 1, is restricted to n < 1. 
However, in contrast to the case of degenerate modes, there 
is no restriction on Gg, as long as U > Gg. For small Gg, the 
entanglement is almost insensitive to U. On the other hand, 
for large Gt , the entanglement depends crucially on the extend 
to which Gg differs from U. When Gg ps U, the correlation 
|?7(b,0)| 2 ~ 1 indicating no entanglement between the modes, 
but | 2 > 1 in the limit of U > Gg. 

These analysis are illustrated in Fig. 4, which shows the 
logarithmic negativity for the case of the two colour entan- 
glement calculated from the previously derived expressions, 
Eqs. (52)-(58), with the drift matrix A now given by 



A = 



( 




V ~1 G (> 





-7 

\Gg 







\Gg 
-7 

-u 



\Gg \ 


u 

-7 



(73) 




FIG. 4. 

forGs 



Variation of the logarithmic negativity Em with n and {7/7 



It is seen that, as before for the degenerate case, an en- 
tanglement occurs for n < 1. However, in contrast to the 
degenerate case, whenever the Cauchy-Schwartz inequality is 
violated, the entanglement occurs over the entire range of U. 



B. Coupled three nondegenerate modes 

We now turn to the problem of determining conditions for 
correlations and entanglement between the three modes. We 
are particularly interested in the problem of entangling two 
independent modes that are simultaneously coupled to a third 
intermediate mode. This situation is encountered in Eq. (69) 
where two mutually independent and nondegenerate in fre- 
quency polaritons are simultaneously coupled to the mechan- 
ical mode. We note that both polaritons are coupled to the 
mechanical mode through a parametric interaction that can 
create an entanglement between the polaritons and the me- 
chanical mode. An interesting question then arises whether 
this kind of the interaction could result in an entanglement 
between the polaritons. We solve Eq. (69) for the steady- 
state and see whether the anomalous cross correlation function 
(6O6H), necessary for entanglement between the polaritons, 
is different from zero. 

To keep the mathematical complications to a minimum, we 
will take <fi + <p = 7r/4so that the effective polaritons 50 
and SH are in the maximally entangled states. In this case, 
Gg = — Gtt = Gt- The situation of <f> + cp = it /A holds for 
the case of the cavity field frequency on resonance with the 
exciton frequency (5 = 0) and the coupling of the exciton 
mode to the cavity mode much stronger than the coupling of 



12 



the mirror to the cavity mode, i.e. ^> G q . 

Let us first examine the stability conditions for the steady- 
state solutions of Eq. (69). It is not difficult to find that the 
eigenvalues of the 3x3 matrix, appearing in Eq. (69), are 



Ai = 7, A 



2,3 



(74) 



where B = ^U 2 - \G\. Since U > G t /\/2, we have that 
the parameter B is a real number for all values of the mechan- 
ical constant Go- It is clear that the transformed operators are 
damped with the rate 7 independent of Go- In other words, 
there is no threshold for Go, which means that the three-mode 
system will decay to a stabile steady-state independent of Go- 
If the mirror oscillates in a thermal field we then obtain 
from Eq. (69) the following nonzero steady-state correlation 
functions 



(SbUb) = n 

(seUe) = (ot^ii) 

(SQSb) = (6115b)* 



(n + 1) [3G t 2 + 8(B 2 + 7 2 )] G 



D 



(n + l)[3G 2 + 8(B 2 + 7 2 )]G 



2D 



i(n+l)(7 + i[/)[3G 2 +4(S 2 -2 7 2 ) + 12i7C/]G t 



D 

{5 ^ gn) = (n + l)[3G 2 +4QB 2 -2 7 2 )-12 t7 ;7]G 



(75) 



with D = 8(B 2 + 7 2 )( J B 2 + 4 7 2 ). 

Unfortunately, the anomalous cross correlation function 
(<50(5II) that is necessary for entanglement between the po- 
laritons, is equal to zero. Hence, the polaritons remain uncor- 
rected and therefore cannot be entangled. In that case, the 
indirect two-photon coupling between the polaritons which is 
provided by the oscillating mirror is effectively zero. The rea- 
son for it is that the oscillating mirror couples to the polari- 
tons with the opposite phases, as it is evident from Eq. (69). 
In physical terms, the oscillating mirror establishes the phase 
difference <f>g — (p^ rather than the phase sum 4>@ + <f> v be- 
tween the two polaritons. In this case, that mutual behaviour 
of the polaritons and the mirror does not create two-photon 
correlations between the polaritons. 

In order to gain some appreciation of the magnitude of the 
first-order coherence, we evaluate the degree of the first-order 
coherence between the polaritons, we find 



|7(e,7r)| 



|3G 2 +4(5 2 - 27 2 ) 4-12i 7 J7| 
3G 2 +8(B 2 +7 2 ) 



(76) 



Thus we see that the oscillating mirror, although not being 
able to entangle the polaritons, it in turn makes the polaritons 
partly coherent. Recall that the polariton modes are in ther- 
mal states, g^ = ' = 2. Thus, two independent thermal 
modes can be made by the oscillating mirror mutually coher- 
ent and the degree of coherence can, in principle, be as large 
as unity. Note that | 7 (6/ j7r )l is independent of n. The reason 
for this is that the coherence depends on the phase relation be- 
tween the modes and the information about the phase is not 



carried out by the thermal fluctuations. Moreover, the expres- 
sion (76) shows that the polaritons are not perfectly coherent, 
the coherence raises with the coupling strength Gt and be- 
comes unity only in the limit of Gt ^> Cl, 7. 

Since the anomalous cross correlations were not generated 
when the mirror was oscillating in the thermal field, we now 
turn to evaluate the cross correlation function (50(511) assum- 
ing that the mirror oscillates in a squeezed vacuum. We have 
already seen that in the two-mode case, the coupling the mir- 
ror to a squeezed field resulted in the generation of the corre- 
lation functions that were zero in the thermal field. Following 
this observation, we recalculate (50511) using Eq. (61), and 
find that in the squeezed field the anomalous cross correlation 
function is now different from zero and is of the form 



[3G 2 + 8(B 2 +7 2 )] . (77) 



The correlation functions determining the occupation of the 
polariton modes (SQ^SQ) and (JIFJII) are not changed by 
the squeezed vacuum and are identical with the results found 
for the thermal field, Eqs. (75). This leads to a quite simple 
expression for the degree of the anomalous cross correlations 



\V(e,^)\ 



1 



(78) 



This expression indicates that the polaritons could be entan- 
gled as long as \m\ 7^ 0. We note that |?j(0 i7r )|, which is 
necessary for the generation of entanglement between the po- 
laritons, is independent of the cavity and the mirror parame- 
ters. It is solely determined by the parameters of the squeezed 
vacuum field in which the mirror oscillates. The reason for 
it is in the broadband nature of the squeezed vacuum, we 
have assumed here, that all modes of the vacuum are equally 
squeezed independent of frequency. 

It should be stressed that the squeezed vacuum creates not 
only the anomalous cross correlation between the modes but 
also creates the anomalous autocorrelations in the modes. A 
straightforward calculation shows that 

(SO 2 ) = (5H 2 ) = ^ [3G 2 +4(B 2 -27 2 ) + 12 l7 ^ , 

(79) 

and then the degree of anomalous autocorrelations are 



1*7(0,0)1 



|»?(7r,7r)| 



|3G 2 



4(£ 2 - 27 2 ) + 12i7C/| 



3G 2 + 8(£ 2 4-7 2 ) 



(80) 



Hence, the creation of the anomalous cross correlations be- 
tween the polaritons is accompanied by the creation of the 
anomalous autocorrelations in the modes. It is easily verified 
that the inequality l^e,^ | 2 > |?7(e,e) 1 1^,^)1 = l'7(e,e)| 2 is 
always satisfied regardless of the value of Gt- 

Normally, we would expect that the requirement for 
the anomalous cross correlations to satisfy the inequality 
IVid.Tr)] 2 > l 7 7(e.e)lh(7r.7r)l is sufficient for entanglement be- 
tween the modes. This was suggested by the inequalities (51) 
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and (72) which were sufficient conditions for entanglement 
between one of the polaritons and the mechanical mode. How- 
ever, substituting Eqs. (76), (78) and (80) into Eq. (41), we 
find that the Cauchy-Schwartz parameter is of the form 



find the following nonzero correlation functions 



X(«,w) 



1 + 



|m| 2 
(n+1)' 



(81) 



where 



ITT = 1 - 



24(G 2 



2B 2 )(B 2 +7 2 ) 



[ZG 2 + %(B 2 +1 2 )X 



(82) 



Since 



1 = n(n + 1) < (n + l) 2 and |T| 2 < 1, the 
Cauchy-Schwartz inequality cannot be violated. We therefore 
conclude that the presence of a strong first-order coherence 
between the polariton modes prevents the modes from being 
entangled. 

We may summarize that the parametric interaction between 
the polaritons and the mechanical mode rules out the creation 
of entanglement between the polaritons. The effect of these si- 
multaneous parametric interactions is to create the first-order 
coherence between the polaritons. 

Consider now a different scenario. Suppose that one of the 
polaritons is coupled to the mirror by a parametric interaction, 
but the other one is coupled by a linear-mixing interaction. 
We shall demonstrate that this kind of coupling can occur in 
our system and it turns out is sufficient to create entanglement 
between two modes that are not directly coupled to each other. 
For example, if instead of the two polaritons, we consider their 
linear superpositions 

8A 1 = —(SQ- 8IL) , 8A 2 = -f= (SQ + SU) , (83) 
v 2 V 2 

we readily find that Eqs. (69) can be transformed into equa- 
tions 

8& = _ 7<5 6t + -^GtSAx + v^yf (t), 
V2 



1 8A 1 - iUSA 2 - -=G t 5V + y/2jA?(t), 
V2 



<L4i 

SA 2 = -7<L4 2 - iU8A x 



(8bHb) 
(8A x 8b) 
(SA 2 5b) 
(SAlSA,) 
(5A{SA 2 ) 
{SA\SA 2 ) 



(n+1) [3G 2 + 8(i? 2 + 7 2 )] G 



D 

-2t 7 (n + 1) [3G 2 + 8(B 2 + 7 2 )] G t 
D 

V2{n + 1) [3G 2 + 4(B 2 + 7 2 )] UG t 
D 

(n+l)(i? 2 + 4 7 2 -lG 2 ) G 2 
D 

6i (n + l) 7 C/G 2 



D 

6(n + l)U 2 G 2 
D 



(85) 



It is easy to see that the anomalous cross correlation func- 
tion (8A 2 8b) is nonzero indicating that the combination of 
the parametric and linear-mixing interactions between three 
modes may result in entanglement between two modes that 
are not directly coupled to each other. Because the anoma- 
lous autocorrelations and the first-order coherence are zero, 
\V(A 2 ,A2)\ 2 = \V{b,b)\ 2 = l7(A 2 ,b)| 2 = 0, the sufficient condi- 
tion for entanglement between the modes 8b and 8A 2 is that 



\V(A 2 .b)\ 2 satisfies the inequality \qi 



A 2 ,b)\ 



> 1. 




2 1 Af{t). 



(84) 



It follows that in this case not 8b^ but 8A\ plays the role 
of the intermediate mode, the mechanical mode 8tf is cou- 
pled to 8A\ by the parametric amplification process, and the 
mode 8A 2 is coupled to 8A\ by the linear mixing process. 
There is no direct coupling between the modes 8b^ and SA 2 . 
Nevertheless, we shall demonstrate that the modes Stf and 
8A 2 can be entangled. In order to illustrate the idea, we pro- 
ceed to determine the anomalous cross correlation function 
(8A 2 8b) from the equations of motion (84). With the help of 
Eqs. (21) and (22), we solve Eqs. (84) for the steady-state and 



FIG. 5. Variation of the logarithmic negativity En with n and U 
illustrating the entanglement creation between the indirectly coupled 
modes 8A2 and 8b for Gt = 7. 

Figure 5 shows the logarithmic negativity En for fixed 
Gt = 7 and gradually increasing n and U. As in the case of 
two coupled modes, the logarithmic negativity is determined 
from the steady-state correlation matrix, which in the present 
case is a 6 x 6 matrix derived by solving an equation that is 
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formally identical with Eq. (55), but with the drift matrix A 
given by 



(86) 



and with the diffusion matrix D = diag[(2n + 1)7, (2n + 

1)7,7,7,7,7]- 

It can be seen from Fig. 5 that entanglement occurs for 
n < 0.5 where it is present for all U. The degree of the en- 
tanglement increases with U, the strength of the linear mixed 
process. Note that the entanglement is reduced when U is 
comparable to G t , and increases as U departures from Gt- 

We conclude this section that two independent modes may 
become entangled by a suitable coupling to an intermediate 
mode. The entanglement occurs when one of the modes cou- 
ples to the intermediate mode by the parametric interaction 
process whereas the other couples by the linear-mixing inter- 
action. 



VI. EXPERIMENTAL CONSIDERATIONS 



Finally, we examine parameter ranges in which the pre- 
dicted coherence and entanglement effects could be observed 
with the current experiments. As a check on the validity of 
the limitation to only the k = 1 exciton, we evaluate the for- 
mula (13) for fk using experimentally realistic parameters of 
an optical lattice composed of 85 Rb atoms [27, 28]. By taking 
N = 10 3 sites, the cavity mode volume V = 10~ 10 [m 3 ], 
the atomic transition dipole moment /i = 5x 10~ 29 [Cm] of 
the hyperfine transition 5 2 Si/2— 5 2 P3/2 in an 85 Rb atom, 
and the cavity frequency u) c on resonance with the first k = 1 
exciton mode that is comparable to the atomic transition fre- 
quency of uj a = 2.5 x 10 15 [Hz], we obtain for the coupling 
strength of the first k — 1 exciton, fi/h= 1.6xl0 8 [Hz], and 
for the third k = 3 exciton, fz/h= 5.3 x 10 7 [Hz] that is one 
order smaller than the k = 1 coupling strength. This simple 
estimation of the coupling strength shows that fk decays very 
fast with k. Thus, our theory can describe the single exciton 
system quite accurately. 



The stability condition for the three mode coupling case 
requires U > G t /V2, where U = J& + G\ and G t ~ Go- 
Since tl ~ |/i \/h = 1.6 x 10 s [Hz] and the typical coupling 
strengths G /2ir w 10 6 [Hz], the inequality U > G t /y/2 can 
be satisfied with the realistic experimental parameters. 



VII. SUMMARY 



In summary, we have presented an analytical study of co- 
herence and correlation effects produced in a single-mode 
nano-mechanical cavity containing an optical lattice of reg- 
ularly trapped atoms. The system considered is equivalent 
to a three-mode system composed of two polariton modes 
and one mechanical mode. We have shown that the system 
is capable of generating a wide class of coherence and cor- 
relation effects, ranging from the first-order coherence, the 
anomalous autocorrelations and anomalous cross correlations 
between the modes. We have been particularly interested in 
the relationship between the generation of entanglement and 
the first-order coherence in the system. The results show that 
the generation of the first-order coherence between two modes 
of the system is equally effective in destroying entanglement 
between these modes. There is no entanglement between the 
independent polariton modes when both modes are simulta- 
neously coupled to the mechanical mode by the parametric 
(squeezing-type) interaction. There is no entanglement be- 
tween the polaritons even if the oscillating mirror is damped 
by a squeezed vacuum field. The intermediate mechanical 
mode effectively creates the first-order coherence between the 
modes. Finally, we have shown that in order to effectively en- 
tangle independent modes, in this system, one of the modes 
should be coupled to the intermediate mode by a parametric 
interaction but the other mode should be coupled by the linear- 
mixing (beamsplitter-type) interaction. 
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